PROTON-NEUTRON OCCUPANCIES OF SPHERICAL 
SHELLS IN THE GROUND STATES OF SOME DOUBLY 
EVEN 2S-1D SHELL NUCLEI 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 

MASTER OF TECHNOLOGY 



f »,j, 

\d 


by 

MALAY KUMAR NANDY 


to the 

nuclear engineering and technology programme 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

JANUARY, 1987 



-? etc 1997 



NBTP " /9S7-A] _ f\jp|\/__ pf^Q 



CERTIFICATE 


This is to certify that the work, entitled 
' Proton-Neutron Occupancies of Spherical Shells in the 
Ground States of some Doubly Even Nuclei 1 has been carried 
out by Mr. Malay Kumar Nandy under my supervision. 

To the best of my knowledge this work is 
original and has not formed the basis for the award of 

wMte 

any other degree. 


St' — * 

Dr. S.K. Sharma 
Asst. Professor 
Department of Physics 
I.I.T. Kanpur. 



(ii) 


ACKNOWLEDGEMENTS 


I wish to express my sincere gratitude to 
Dr. S.K. Sharma for his constant guidance and encouragement 
at all stages of the work. 

I wish to thank Prof. K. Sri Ram, Prof. A. 
Sengupta, Dr. M. S. Kalra and Mr. P.Munshi for their help 
and useful suggestion. 

I thank Dr. G. Mukherjee for his help 
regarding the computational work. 

Finally, I wish to thank Mr. Ajay Krishna Roy, 
Mr. Amit Sircar and all other friends for their constant help 
and encouragement. 


MALAY KUMAR NANDY 






(iii) 



CONTENTS 




PAGE 



CERTIFICATE 

(i) 



ACKNOWLEDGEMENTS 

(ii) 



LIST OF TABLES 

(v) 



LIST OF FIGURES 

(vi) 



ABSTRACT 

(vii ) 

CHAPTER 

1 

INTRODUCTION 

1 



Reference 

4 

CHAPTER 

2 

THEORETICAL FRAMEWORK 



2-1 

The Shell Model and the Residual 

Interaction 

5 


2-2 

The Hartree-Fock Method 

8 


2-2-1 

The choice of an expansion for 

the HF orbits 

8 


2-2-2 

The symmetries of the HF field 

9 


2-2-3 

The HF equations 

11 


2-2-4 

Expansion of the occupied orbits 
in the HF equations 

14 


2-2-5 

Single major shell HF calculations 

15 


2-2-6 

The use of a reference closed shell 

nucleus 

16 


2-2-7 

The iterative procedure 

17 


2-3 

The Hartree-Fock-Bogoliubov Method 

19 


2-4 

The Technique of Angular Momentum 
Projection 

26 



(iv) 


2-5 Calculation of Subshell Occupation 
Numbers in Terms of the State with 
Good Angular momentum Projected from 
the HFB Intrinsic States (The Projected 


Hartroe-Fock-Sogoliubov Method) 29 

2-6 Effective Interactions 32 

2-7 Calculational Procedure 34 

References 35 

CHAPTER 3 RESULTS AND DISCUSSION 


3-1 A Comparative Study of the PHFB and 
the Shell Model Approaches 
Employing the Nildenthal Interaction 36 

3-2 A Comparative Study of the PHFB 

Predictions for Proton and Neutron 
Subshell Occupation Numbers 
Resulting from Various Effective 
Interactions Operating in the 
2s-ld Shell 37 


Ref erencos 


51 



TABLE 

TABLE 

TABLE 


LIST OF TABLES 


I Ground State Subshell Occupation 

Numbers from Experimental Data 

i'T The Subshell Occupation Numbers 

for Protons and Neutrons Resulting 
from Various Effective Interactions 

III Comparison of the Subshell Occupation 
Numbers Resulting from the PHFB as 
well as the Spectral iVethod fr c 
Some N=Z Nuclei in 2s-ld Shei. 


(v) 


PAGE 


39 


43 


46 



LIST OF FIGURES 


PAGE 


FIG. 1 

FIG. 2 

FIG. 3 


FIG. 4 


T-l Matrix Elements of Various 
Effective Interactions 47 

T=0 Matrix Elements of various 

Effective Interactions 48 

Comparison of the PHFB as well as the 
Shell Model Predictions Resulting 
from the Wildenthal Interaction for 
Some 2s~ld Nuclei ^ 

Comparison of the PHFB Predictions 
Resulting from Various Effective 
Interactions Operating in the 2s-ld 
Shell 50 



(vii) 


A BSTRACT 


The aim of the project discussed in the thesis 
is twofold. First we assess the relative efficacy of the 
Projected Hartree-Fock-Bogoliubov (PHFB) approach for 
nuclear many-body problem. This is done by comparing the 
subshell occupation numbers resulting from the PHFB as 
well as the exact (configuration mixing) shell model 
approaches. Second, we discuss the calculation of subshell 
occupation numbers for a number of 2s-ld shell nuclei employ' 
ing various effective interactions in conjunction with the 
PHFB method. The task is facilitated by the fact that the 
PHFB approach is considerably faster than the usual (shell 
model) configuration mixing method. A comparison of the 
calculated results with the available experimental infor- 
mation yields significant information regarding the suita- 
bility of various interactions for detailed microscopic 
studies in the 2s-ld shell. 



CHAPTER 1 


INTRODUCTION 

The expectation value of the subshell occupancy- 
operator is one of the simplest measurable dynamical 
quantities in the ground state of a nucleus. The measured 
values of the subshell occupation numbers provide a direct 
test of the nuclear wave functions involved in the micros- 
copic description of the ground state. 

A calculation of the subshell occupation numbers 
in terms of the exact shell model wave functions is not 
feasible for the 2s-ld shell nuclei; the dimensionalities 
of the Hamiltonian matrices involved are prohibitively 
large. Furthermore, the collective nature of the observed 
properties of the nuclei with 16<A<40 necessitates the 
involvement of the full 2s-ld valence space, and the 
truncation of the shell model space is not expected to 
give accurate results. 

Potbhare and Pandya^ sometime ago reported the 

calculation of the ground state energies and the subshell 

occupation numbers of several 2s-ld shell nuclei in the 

framework of the spectral distribution method developed by 

2 3 

French and co-workers. 5 The occupation numbers resulting 
from a number of effective interactions in the 2s-ld shell 
were compared. However, it turns out that the ground states 
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in various 2s-ld shell nuclei are usually more than five 
widths (o) away from the centroid of the energy distri- 
butions. This severely constrains the usefulness of the 
spectral distribution method involving just the first two 
moments of the Hamiltonian; the results obtained via the 
spectral methods become very sensitive towards the 
higher-order shape correction to the configuration level 
density. 

In contrast to the spectral methods, the Hartree- 
Fock (HF) or the Hartree-Fock-Sogoliubov (HFB) approaches*^ 
are a priori more suited for the calculation of subshell 
occupation numbers for the ground state resulting from a 
given effective interaction. The requirement of the 
minimization of energy for the HF or HF3 ansatz of the 
variational wave function is connected primarily with the 
extraction of the one-body field. The lowest-energy 
HF/HFB intrinsic state, therefore, corresponds to the opti- 
mum single-particle field (resulting from two-body intera- 
ctions) appropriate for the ground state. This single- 
particle field in turn provides the best set of (single- 
particle) wave functions for calculating the expectation 
values of various one-body operators in the ground states. 

The aim of this paper is twofold. First we compare 
the subshell occupation numbers resulting from the projected 
Hartrec-Fock-Bogoliubov (PHFB) and the shell-model approach 
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employing latest 2s-ld shell interactions; the comparison 
permits an assessment of the relative efficacy of the 
PHFB approach. Second, we report on the calculation of 
subshell occupation numbers for a number of 2s-ld shell 
nuclei employing the angular momentum projected HFB wave 
functions resulting from various effective interactions. A 
comparison of the calculated results with the available 
experimental ones provides information regarding the suita- 
bility of the effective interaction vis-a-vis a detailed 
microscopic study of various isotopes. 

In Chapter 2 we give an outline of the formalism 
for the calculation of the proton and neutron subshell 
occupation numbers using variational HFB wave functions. 
Chapter 3 contains discussions on the theoretical results 
and their comparison with the available data obtained 
from various pickup and stripping reactions. 
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CHAPTER 2 


THEORETI CAL IR1AMEWGRK 


The Shell Model and the Residual Interaction 


Normally a nucleus is a many body system, so we 
are concerned wi th solving the many body Schrodinger equation 

H ( 4 ) = 3 (jj ( 1 ) 

where (j) is the intrinsic state of the nucleous assume!"' to be 
given by the Slater determinant: 

^ ^;j 1 3 2 ...,G A (x 1 ,x 2 ,.*.x A ) 

= A- y ( X l) (*9) **• X A ) 

vaT ^ ' 1 " 

^(3 9 ^ X 2^ *** ^p 2 ^ x A^ 


^,/ x l^ C ^p A ^ x 2- ) *** ^ a ^ x A^ I 

(2) 

(j),, (x.) being the single-nucleon wave functions of the orbit 
Pi 3 

The coordinates Xj include description of position, spin 
and other coordinates the wave functions mav be dependent on. 

The many-body nuclear Hamiltonian in eq.(l) is given 


+ * V ij 

i<J J 
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where A = Mo. of nucleons present in the nucleus 
r •= K. £. of the i 1 nucleon 
V^i= V(x^,Xj) = The a ctu al two nucleon interaction 
potential 

The two-nucleon interaction is vry strong and 
contains a short-range strongly repulsive core in some two- 
body states. .So, its form is too much complicated and yet 
the exa ct form is unkno wn . Thus it is hardly possible to 
have an exact solution of eq. (1). 

The success of the independent particle model in 
explaining a large number of experimental data suggests 
that tho short-range rupulsive core does not play an 
important rolo in many nuclear properties. It is therefore 
believed that the wave functions which are the solutions of 
eq. (l) with an effe ct ive nucleon-nucleon interaction from 
which tho repulsive core has been removed are significant. 

In an attempt to solve the many body eq. (1) 
approximately, we assume that each nucloon of the nucleus 
.moves independently in a potential field that is average 
of the interactions of this nucleon with all the other 
nucleons of the nucleus. Denoting this one-body interaction 
by v i = v(x i ) we have 

A A 

H = E t- + Z 
i=l i=l 


v . 
1 


( 4 ) 
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dut in writing so vie are neglecting a 

when (3) is written as 

A A A 

H = L t • + l v • + ( 1, 
i i=l i=l i <j 


term, 



which is obvious 


2 v.) (5) 

i=l x 


be write the quantity in parentheses as 


A 

( £ 

i<j 



A 

E 

i=l 


Vi) 



res 



res 

where 'A ^ is the residua^ two-nucleon interaction that 
remains after the extraction of the one-body interaction 
from the effective two-nucleon interaction. 


The correctness of the mathematical results of (4) 
depends upon how we choose the form of v... 

In shell-model calculations we choose the form of 

v i as 

v i = I M L ' 2 f l + K + T Vh (6) 

and set =0 for all i , j . The oscillator parameter (to) is 

chosen to match the observed r.m.s. radius of the nucleus 
and the strength of the spin-orbit coupling (a,y) are chosen 
to match the spin-orbit splitting. 

The choice (5) permits one to mate a number of 
useful prediction for core + one-particle systems. For 
nuclei involving a few valence particles one can attempt to 
diagonalize V ^ exactly in the basis provided by (6) . 
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However, for larger (n>3) number of valence particles this 
procedure becomes quite cumbersome; the dimensionalities 
of the rolcvent Hamiltonian matrices become prohibitively 
large. In such cases one invokes tne variational methods 
such as the Hi- or the Hi-3 methods. 

2-2 The Hartree-Fock I.:e tho d ^ ^ 

In the Hartfee-Fock formalism we denote the 
correct single-nucleon states } p> be expanded in terms of 
the basis states Ji> 

jp = I 1 i > ; 1 < p < A (7) 

i 

The basis states may be taken to be the harmonic oscillator 
states ( in tine presence of spin-orbit coupling ) used in 
ordinary shell-model calculations. 

o-2-l The Choice of an Expansion for the HF-orbits |j3> 

The basis |i> on which HF-orbits | h '> are expanded 
is arbitrary in principle. In practice, useful expansions 
are suggested by 

i) the convergence of the expansion, 

ii) the convenience of calculating the matrix 

elements of the two-body interaction, and 

iii) the symmetries of the HF- field 

Harmonic oscillator orbits were found to give 
convergence 

excellent/for the expansion of bound orbits of spherical 
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closed -shell nuclei. Thus harmonic-oscillator wave functions 
can be used for the radial part of the basis state ji>. 

Shell-model states }nljm'f> are convenient to use 
because of the general availability of matrix elements of 
the two- body interaction in this representation. 

Expansion (7) for the Hr-orbits on the shell-model 

basis is 

IP =nl? m ^ |nlj "V (8) 

Summation over T is avoided to ensure that proton and 
neutron states do not mix together; T is the third 
component of iso-spin t, distinguishes proton and neutron 
states . 

2-2-2 Symmetries of the HF-Field 
Axial Symmetry: 

The HF-orbits are said to be axially symmetric if 
each orbit | j 3 2 has a good projection of r.e. one orbits 
are eigen states of J z 

J z iP> = If> 

Thus each state 1 p> is characterised by a given projection 

m 0 and hence the expansion of an orbit p is limited to 
P 

basis states of a given 

|P> = C nl j | nl j m p r p> 

nig J 1 


( 9 ) 
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Obviously such orbits would correspond to a potential 
which has rotational symmetry about z-axis. Also the total 
wave function (j' given by (2) will have a definite projection K 

J dV = K iV 

where K = m n +m D + .... + m_, 

d f>A 

Time-Reversal Symmetry: 

With real . coefficients, time reversal 
symmetry is equivalent to the symmetry of a rotation of 
71 about the y - axis ( equivalent to a rotation in the xz-plane) 
times the parity operation 

T = e y P 

In even-even nuclei, time-reversal symmetry is 
obtained by requiring that, for each occupied orbit 

= nljm 

the time-reversed orbit 

l?> =T|g> = S <^ 11n (-D 1+j " m inlj,-mT 6 > 

is also occupied and orthogonal to j|3>« Then the HF-field 
commutes with the time-reversal operator and the orbits |p> 
and i p"> are degenerate. Also 

T | (|)> = \$> 
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Obviously, in axially symmetric solutions the 
orbits jj3> and ( |3> are always orthogonal. In other 
solutions they are orthogonal if all the components of the 
orbits have either 


m+l+l/2 even 

or m+1+1/2 odd 

i.e. expansion of each orbit is limited to basis states with 
either m+l+l/2 even or m+l+l/2 odd. 

2-2-3 The Hartree-Fock Equations 


Wo 


consider axial -symmetry so that each orbit 


p> has a definite projection rru and is expanded as 

|,C> = 2 4 

3 p 


( 10 ) 


where j stands for nlj and m. for m,T-. This expansion is 

P 1-' P 

subject to the normalization requirement 


2 |C? j 2 = 1; for all p 

j 

The total nuclear state, given by the determinant in eq. (2) 
then can be labelled by the total projection K 

i. 0 . j(jv> = YXf det i Pg > Ip a >3 


The energy of the system is obtained by taking 
the expectation value of the Hamiltonian given by (3) in 
the state |(j)^> t 

E hf = «t K |H|(t) K > 

A A 

= 2 o|t|p> + 2 <j3jU|Vjpp.> 

p=l 
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wno re i-ptp is called the HF-energy of the system, <j?jtjp> 
is the U.fc. of a nucleon in the orbit j and <pp,j V| j3p> is an 

an ci symmetrised matrix element of the two-body interaction 

v . . . 

1 j 


On expanding the orbits |j5> as given by (10) we 


have 


■'-LTL' 

iii 


= < 

* 
3 1 

= 12 C, X 

P jj.J 3 P 


C.f 

P 


<j 1 m e |t|j 3 m,> 


+ 2 
P<U 




h> 


( 11 ) 


flow, the purpose of the HF theory is to determine 
the single-particle orbits |p> (eg. 10) in the Slater determi- 
nant of This is equivalent to the problem of determining 

the C] coefficients in eg (10), since |jm..> are known basis 
P P 

states (chosen to be shell-model states). 


This is achieved by requiring the energy of the 

system ( i.e. 2j_jp) stationary and minimum for infinitesimal 

variations of the urbits p ( brought about by infinitesimal 

variations of tho C-2 coefficients) subject to the constraint 

P 

given by eq (10). Thus we require 




2 

h 



c p 1] 


0 


( 12 ) 


[«y H i<k> - 
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where €p arc introduced as Lagrange multipliers. V/ith the 
help or the expansion given in eq (11), this derivative can 
bo calculated directly: 


^1 3 j 0 

t V I j 3 m b > + 2 Z <j,nu; Mvjjmiu; \i> 

3 3 ' ■ ' \i j 3 p x ° ' L 

j'l 

— £. C - 1 

- V p 

which is the result after removing the common summations over 
p. This can also be written as: 


T r |t| j q m p > + 2 pjvjjmiu; p>] C 

4 ~ L W ° l- 1 , n — <j O u 


^3 


|a=l 


= C i 


'1 


(13) 


which has the form of an oigen value equation 

3 3 „ ii 


2 < j x mp | h | j 3 ny C, = * 
J 3‘- 


(14) 


whore 


A t 

<jl0.j|h| j 3 m ]? > = [ <Jj _ m 3 1 1 1 j 3 m s > + 2 <0^; P- { V | j 3 m^ ^ 3 
1 1 ' h 1 


(15 ) 


h is called the HF hamiltonian for a single nucleon 
and the ouantities are tho HF single-particle energies. 

' r‘ 

The states jp> are the HF single particle states and can 
be known on determining the coefficients. 

The equations (14) and (15) are the HF equations 
for the orbit j3. They are solved by iteration procedure. 
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2 2~4 Expansion for the Occupied Orbits in the HF Equations 


Under axial symmetry, the p-orbits in eq (15) can 
be expanded as 


V<> 


* c l hv 


whence the two-body term in eq. (15) becomes 

^ 2 r ?4 / . 

V S <J i*V J 2 m p* v 1 


.* 

I ? S c u 2 C <Jx»«* V«l v lj.w PV 
J 2 3 4 


+ \ l ( ijy C '/ }< W t m u > 

Thus we can not separate out the factor 

P = r c^ 2 c 4 
l h_U n » ** 

A 

and hence can not define the density matrix p « However, if 

t 

we violate axial symmet ry in the expansion of |p> i.e. 

jm 

!p> = 2 C |jm> 


jm 


u 


we can define density matrix as shown below: 


2 <j,m p ; p | V j j ofii - ; p> 
p 1 ! J J y 


— V 

— o 


2 ( 2 C 


j 2 m * c J 4 m ' ) < jjmQ f j 2 m { V | j 3 m 0 ; j 4 m 1 > 


j 2 m j 4 m* i 


, P P 


(16) 


as we can separate out the factor in parantheses. They 


are matrix elements of the density operator 
p A-16 j„m* j 4 m' 

= - c ,, c 


r- 


p^=i 


where jm stands for nljmT. 
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The uncoupled, two-body interaction matrix elements 
in (lo) are given in terms of anti symmetrized, coupled two- 
body .matrix elements by 


<j ituT 0 ; j 9 mljv| j 0 m,- T-. ; 

J- Jw 1 «— O .J , 


= t v r 1 ^ i r 

tt r fn ■ *11 -> m is. llln 

1\. X i '‘irj-. p P 


^3 ^ 4 


J 1/2 1/2 I 1/2 


m’ K ^ T T i'.ijt 7. 


1/2 T -I 

T> ?, 

i -“T 


x <(j 1 j 2 )Jljv! (j 3 j 4 )JT> 


(17) 


The quantities in square brackets • are the Clebsch-Gordon 
coefficients s 


[ 1 2 ] = < J’2 m 2i ^1 J 2^ JM> 


J 


rrij m 9 M 


vdaich 


are non zero only when m^+nip=iv; and hence suggest that 


rn= m' 
T= T' 


2-2-5 Single iuajor Shell Hartree-Fock Calculations 


In single major shell HF calculations the expan- 
sion (8) of each HF orbit is limited to states | nljmT> 
belonging to a major shell of the spherical oscillator. Thus 
for (ld^yp 2Sjyg ld 3/2^ ma J or shell calculation we expand the 
HF orbits as 


! iB> = 


nljrn^(ld,2s) 


'nljm 


nl jmT„> 

P 
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for the valence nucleons in the (id, 2s) major shell. The 
inner shells which are completely filled are assumed, to 
form snjinert .jsojrsj , and are ignored, in actual calculations. 
Also A should be replaced by (A~I!o. of nucleons in the core) 
= ( A-1S) . 


2-2-6 Tho use of a Reference Closed-Shell Nucleus 

, 16- 

For 2s-ld shell nuclei we assume an inert 0 
closed-shell core plus A— 16 valance nucleons occupying the 
orbits i?. which are expanded on the ! < % / /’2 , ^ s l/2 , ‘ 1 ' <: *3/2 s ^ a ^ es * 
Thus it is convenient to choose ^0 as a reference nucleus. 
Then the sum over \i in eq . (15) can be splitted as: 


rV"" J- 


+ 


£ 


ji=l ii=ls,lp 

'where the first summation is over 
and the twelve 2P— shell orbits of 


<jrV 




p^=l 

therefore Is— shell orbits 
lo 0. eq. (15) then ^becomes 


£ <j.au*n|v| 
ds,lp 1 p 


A-16 

+ £ < J-, m . ; \x V j 3 m P :p.> 

The quantity in parantheses is the Hr Hamiltonian of tne 
closed shell 16 C. Thus the HF Hamiltonian takes the simple 


form 


<j 1 m p |h|j 3 m p 


e. d. 

Jl J l-o 




+ 


£ ^3i m B ' P I M i 3 m p * P^ 
p=l p F 


(18) 
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where the K.E. term in (15) is replaced by the single- 
particle energies e ^ of the ^°C field. These single-particle 
energies may be obtained from the binding-energy differences. 

e j = f 17 °]j - f 16 °] ; je.(2s s ld) 

where [^0] is the binding-energy of ^0, and [^ 7 G]j is 
the binding energy of a state of 17 0 with spin j. Several 
other methods are there for the maipulation of -values. 

An important noint is to note that, the e^.-values 
in (18) are the ones obtained with the °C value of a— '{Ma/h , 
the oscillator constant, As cc varies from one nucleus to 
another, the e^'s also vary. This effect may be taken into 
account by taking parametrized values of e^'s which are 
available. 


2-2-7 The Iterative Procedure 


Finally we have the HF equations, fox' the 2s~lc! 


major shell calculations: 


1. Z | hj j 3 m 0 >CC ° = C/ (19) 

j 3 ‘ 1 p 

A— 16 

2. |h| = 6^ dj ij3 + ^ <j 1 m p ^!v|j 3 m ? j (O 


A- 16 j 9 m* j.m 
= e. d. . + Z E ( -Z C C ) 

J i J 1 J 3 j ? m j 4 m F— 1 ^ 


<J l ffi p ; ^2 m i V l j 3 r V J *4 m> 


( 20 ) 
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for the orbits p. They can be solved by iteration as 
follows : 

a. An initial guess for the set of coefficients is 

h 

made . 

b. V/ith this set of coefficients , the matrix elements 

of the HF Hamiltonian are calculated from eq (20)* 

c. Knowing the- HF Hamiltonian matrix in the step b. , it is 

diaaonalised to solve the eigen-value eo (19). € 1 s and a 

* p 

now set of coefficients C?. m are obtained. 

Now we return to step b. This process is repeated until 

successive di agonalizations produce an almost equal 

sot of coefficients C^ m ( within desired degree of 

P 

accuracy) . 
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2—3 The Hartree-Fock-Bogoliubov 


i lethod 1 


VJe consider the HF field possessing time-reversal 
symmetry. For even-even nuclei this symmetry is easily 
ensured by requiring that for each occupied orbit |k>, the 
time-reversed orbit |k> is also occupied. 

The axially symmetric intrinsic deformed HFB state 
with K=0 can be written, according to the BCS theory of 
superconductivity, as 


IW- = l <V V k b k b k ) !°> (2D 

This amounts to putting a pair in the states k and k with 
probability (Here and indicate respectively the 
fullness and emptiness of the orbit k). Thus |(j) > contains 
a superposition of many pair-excited configurations each 
having different number of nucleons. In the case of the 
axially symmetric case, the creation operators can be 
expressed as 


? C k a jk 


4 = ] c k a l-k 


( 22 ) 


aT^ creates a particle in the state |jk>< 


The normalization of the state |(j) > requires 


U k + 



1 


( 23 ) 
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In the HF limit V k =1, U k =0 for the occupied orbits 
and V k = 0, U k =l for the unoccupied orbits. Thus 


l fe 

operators b 
represents 


HF 


% 


b, 


k (occupied) 


b± |C> 


(24) 


can perform a canonical transf ormation from the 


k uu a iiisw sot of operators q k such that |(j) > 


to a new 

a vacuum of the new operators q k 


% = u k b k - v k 4 

a -k = u k b ir + v k 

N>o> * f° r a -*-l k. 


(25) 


The Hamiltonian of the many-body system consisting 
of the non-interacting HF particles is given by 

h (HF) = z £ o b + ^ (26) 

k 

where is the HF energy of the level k. 

Though the HF theory extracts the best single-particle 
field from the residual two-body interactions, some parts of 
the latter are still left untreated. Therefore v/e invoke the 
p airing 1 interaction between the particles in the HF scheme; 
this interaction generates the admixtures of pair-excited 
configurations in the intrinsic wave function (j) . V/e write 


the Hamiltoniaa as 
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H = 


S 

k 


£ k K \ 


E 

kk’ 


<V 


V rSS Ikk> bj,b± b^b k 


(27) 


The matrix elements of the pairing interaction 



In what follows we consider the state-independent 
pairing interaction characterized by the average of the 
matrix elements G^,, 


G = <k 'k ' 1 V res i kk> (28) 

We can thus write eq. (27) as 

H = 2 E k b+ b k - G E ^ bj, b±, b k b k (29) 

k kk 

For G=0, this Hamiltonian goes over to the HF limit. As we 
can easily see, for G ^ 0 the pairing interactions smooth 
out the jumps in the occupation numbers for the HF orbits. 
Now the number operator is defined as 

£ = Z b* b k (30) 


In the quasiparticle vacuum 

<d) I n |0 > = n = 2 E V? (31) 

N o 1 1 T o k c 

Since the transformation from particle to quasi*- 
particle formalism does not leave the number operator 



diagonal , we must add to the Hamiltonian (27) a term p n 
and then choose p such that the expectation value of n is 
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equal to the desired number of particles, i.e., we work 
with the Hamiltonian (H-pn). 

Using (29) and (30) we have 

H - pn = I bj b k - G Z bj, b™ f b k b^ (32) 
k kk ' 

where , „o 

= E k - ^ 

i.e. we measure the energies of the particles w.r.t. p (the 
chemical potential ). Later we choose p such that <n> =n. 

After neglecting the terms V k we have 

4 0 |H-pnj(j) o > = 2 S * k V* - G ( S U k V k ) 2 (33) 

k k 

Our task is to minimize (33) : 


d V, «t>o I H “ ‘ un I > 


under the condition 
u k + v k = 1 

Using (33) in (34) we have 


dU k 

4 €,.V k -2G( S U k V k ) (U k+ V k 5 ^) = 0 

On differentiating eq. (35) with respect to V k we get 


au_ k 

5v k - = 



(34) 


(35) 
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therefore we have the set of equation: 

r 2 


V 

2 €, V, - r,( n k 


*k k 


where we have used 


'k u 


■) 


k 


9 o 

+ v k 


= 1 


-s= G 2 U k V k 


In view of eq. (37) the set of eqs. (36), (37) 
solved easily. T /e put 


(36) 

(37) 


(38) 


can be 


V k = cos© 
U k = sin© 

from (36) we have 


or 

or 


2 e k u k v k = A(ug - v 2) 

© k sin2© = -^cos2© 
tan 26 = - (4/ •:, ) 


cos 29 U R V k y-^2 +a 2 


and 

which 


u; 


k 


+ Vf 


together yield. 

TI^ _ .1. ( 1 4- 

u k 2 v 


v k = 


(i- 


\f € k 2 + 4 2 

£k • 

o o 
f 4k~ + 4 a 


(39) 

(40) 


These are called the Hartree-Fock-Bogoliubov equations. 
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At this stage we note that had we started with 
the more general pairing interaction 

G k'k = < k ’k' |v|kk> 

and carried through the minimization of energy we would 
have got more general HF3 equations; 


i [ i + j 


(41) 


I C 1 


f[( E^) 2 +4 2 ] 


(42) 


2 k Y[( E° k ,-,) 2 + A2,] k 


(43) 


also we have 


n = 2 I V( 


(44) 


The iterative procedure for solving the HF3 equations 


We define 


P ( [i, ) ~~ 


f[(F°-p) 2 +^S ] 


(45) 


Q( 


On taking the sum on 


k T[(E°-p) 2 + A? 


(46) 


both sides in eq. (42) and using eq.(44) 


n = N„ 


P(u» + h Q(4 S \) 


we have 
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where N q (=6 
paired states 

V> = 


for the 2s-ld shell) is the number of the 
kk available. Thus we have 
P(H» - N Q +n 




(47) 



1 s <kk|v|k'k'> 
k ’ f[(E k ,“H) 2 +^ k J 



(48) 


The iterative procedure employed in solving the 

HFB equations is as follows: 

a. An intial set of values for p and A^'s is guessed, and 
also which orbits are occupied. 

b. With this set of p and.A k r s eq. (47) is iterated to 
obtain a set of sel-consistent value of p. In this 
iteration theA^s are held fixed at the values 
obtained in the previous step. 

c. y/ith the self-consistent value of p obtained in step b. 
and theA k 's ( which were earlier held fixed) we iterate 
the equation (48) to obtain a self —consistent set of 
values of A k 's. ( in this iteration p remains unchanged). 

d. We return to the step b. 



26 


2-4 The Technique of Angular momentum Projection 


In the axially symmetric HF, each single-particle 
state has a given projection quantum number, and hence the 
HF intrinsic state determinant (b., has a given total projection 

ou an turn number K, which is the sum of individual single- 

* 

particle projection quantum numbers. 


It is usual to describe the actual nuclear state 
with a definite angular momentum J. Therefore the determinant 
(j)^ is written as a linear superposition of states ijijj 




(49) 


Operating with the rotational operator 


ft (it) 


-iaJ I iyJ z 

e e 7 e 


(50) 


where Jl represents the Eulerian angles a,© and y we have 


R(A) |(j) K > ~ 2 3j (Si) 

JM ‘ - 

where are the matrix elements of the rotation operator 

Ms. 

in the basis |^>* 


over XI 


Now multiplying both sides bv and integrating 

where 


2n 

I = / dy 


IX 

f d© sin© 
o 


2tg 

/ da 

o 


3 71 \ -v 

=27+T d jj' °m' 


and using the property 
/ ^0-K (Sl) D M'K (Sl) 
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and finally omitting the primes, we have 

^ / dO D^(JZ) R(J1) |<t> K > (51) 

871 a J 

v/e say, in view of the eq. (51), that the nuclear state of 
angular momentum J and projection M in the Laboratory frame 
(belonging to a rotational band) is projected out of the 
deformed intrinsic HF/HF3 state )(()£>, which has axial 
symmetry and projection K on the z~axis of the intrinsic 
frame. Thus the projection operator is written as 


J 

MK “ 


/ 

aj 


Putting M=K the equation 




R(-Sl) 


I4>k> = 22±i f d n DjJ K (ii) S (Si) i<(i K > 

Stt: aj 

can be represented by 



(52) 


where 

P KK = R (Si) (53) 

K 87t a j 

The projection operator PjL has the following 
properties 


(a) 

II 

P J 

*KK 

(b) 

( P KK> 2 “ 

P J 

*KK 

(c) 

[ H, P KV ] 

= 0 
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The property (a) follows from the fact that angular momentum 
operator is hermitian. Again, the projection operator opera- 
ting on 1^> gives back (l|j^>; therefore property (b) follows. 
Property (c) is obvious as H is a scalar operator and cannot 
change the angular momentum of a state on which it operates. 


In view of the above properties of the projection 
operator and eq (52), the energy of the projected state | l|j^> 

p = < l H l tg > 

reduces to ^ <^jH P^ K j <t> K > ( 54 ) 

Ej ~ «t> K l p jlV 

Using the expression for Pj^ and R (-H-) from (53) and (50) 
respectively and using the fact that exp (-iyJ ) and exp 
(-iaJ 2 ) Operating on the bras and kets <(j)^| and 
produces exp (-iyK) and exp (-iaK) respectively and using 

= e" laK djj 1< (9) e“ iyK 

and cancelling the integrations over a and y we finally 
have 

f n d9 sin9 dj^(©) e |^> 

TC 

/ d9 sin9 d^(9) <(j) K j e~ 10J y |<|) K > 

0 


( 55 ) 
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2“5 Calculation of Subshell Occupation Numbers in Terms of 

States with Good Angular Momentum Projected from the 

HFB Intrinsic States ( The Projected Hartree-Fock- 

4 5 

Bogoliubov Method) 


The axially symmetric intrinsic deformed HFB states 
with K=Q can be written as 


,> = % (U. + V. bL bt™ ) j 0> 


•'O' " . " vw im • im “un im 
im 

where the creation operators bT^ can be expressed as 


(56) 


b im 


? C L a+ jm 


b +- = L cJ m 3 t,. m 

Here, the index i is used to distinguish between different 
stat ' S with the same m, and j labels the spherical single- 

particle orbitals ld 5 y 2 > 2s l/2 ? and ld 3/2* The W3Ve functl0n 
in eq. (56) can be reduced into the form 


IV 


N exp [| Z f ap a* a* 


] 1 0 > 


with 


' a 


8 


V. 

im 


' a{3 


i im a im p U im 


a 


a 


m a ,Hcn p 


(57) 


( 58 ) 


Here a denotes the quantum numbers (j^) and N is the 


normalization constant. 
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Now, the subshell occupation number operator is 


given by 


£> a+ a+ 

d a • 


jm jm 


( 59 ) 


m=- j 


The expectation value of this operator with respect to the 
angular momentum projected HF3 states 

jtiA = p J |(j) > 

1 • 0 00 IY 0 


can be written, in view of the eq. (55), as 

% T 

/ p(9) d~ (9) sin © d© 

TJ i = O 

J n T 

/ n (©) d^ o (©) sin © d© 

where p(9) = «jL | n., e~ 10J y |(j) > 


(60) 


( 61 ) 


n(9) = «J) 0 | e" 


VI V 


(62) 


they arc- related as 


p<8)= n (9 ) [ 2 l° jm , jm ] 

m=- j J y J 

The density matrix ^(©) is defined 

P(q) 

1 ~ T+mI©) 


(63) 


(64) 


whore hi(©) = F(9) f" 


(65) 


and F ap (©) - l d , (g) d P (e) f j m (66) 

m a m {3 a a m p m p “p P 
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AcfRin the normalizations n* aro given by 

syr 

n = / n(6) d^ Q (9) sin 6 d© (67) 

v /h ore n (©) = y"det(l+M(@) j (68) 

It is worth mentioning here that an approximate 
expression for the subshell occupation numbers can easily 
be obtained in terms of the expectation value of the operator 
■fij with respect to the intrinsic state J Cj) o > » 

r) • «b I r) . | (b > =2 I C-? 1 2 V? 

J Y o 1 'j 1 Y o 1 lm 1 lm 

The above equation provides the values for occupation 
numbers averaged over the complete rotational band. These 
band averaged values are likely to differ significantly from 
the estimates based on the use of the projected J=0 wave 


functions . 
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2-6 Effective Interactions: 

In this thesis we shall discuss the results 
obtained with seven alternative two— body interactions. Five 
of these were obtained via least square search so as to 
optimize the agreement between the calculated and the observed 
properties of nuclei in the 2s-ld shell. In what follows we 
give first a short hand level for each Hamiltonian, then the 
number of parameters which were adjusted via least square 
search, and then a brief description concerning the origin 
of the interaction^. 

K+.SPE (no free parameters ) This is the realistic 

effective Hamiltonian derived from the Hamada- 
Johnston potential for nucleon-nucleon scattering 
by Kuo^. 

KB+5PE (no free parameters) This isthe realistic effective 

O 

Hamiltonian derived by Kuo and Brown in the same 
manner as the proceeding interaction. 

K+12FP (9 free parameters) The interaction involves 9 

twc-body matrix elements which were parametrized. 

All other two body matrix elements were held fixed 
at the Kuo values. 

6 

RIP (11 free parameters) The interaction stands for 

11 Radial Integral Parametrization ' 1 . It involves 

11 linearly independent combination of radial 

Talmi integrals which were considered free parameters. 
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M§Dj. (4 free parameters) The letters MSDI^ 1 stand for 

' ' Modified Surface Delta Interaction '* . The 
four parameters are the two surface-delta 
strengths ( one for T=C and one for T=1 inter- 
action), plus two modifying monopole strengths 
( one for T=0 and one for T=l). 


P’:/ 


(18 free parameters) This interaction was 

q 

obtained by Preehom and Wildenthal . Those two- 
body metrix elements of K+SPE which do not involve 
the dg /2 or bit were varied as free parameters. In 
addition, the centroids of the anc * 

®3/o“d3/2 i n teraction were also varied. 


(63 free parameters ) This interaction was obtained 

by Wildenthal^. The 63 two-body matrix elements 

characterizing the full 2s-ld shell interaction 

were obtained by requiring a stable best fit to 

a set of 440 data which sampled every 2s-ld shell 

system with experimentally known energies. An 

overall constraint was the assumption that the 

0 3 

matrix elements decreased with A 3 . s (18/ A) . 

In all the above interactions the single-particle 
energies were also considered as free parameters. 
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2-7 Calculational Procedure 


The initial guess for the wave functions involved 
in HF/HF3 iterations was generated bv diagonalising the 
Nilsson Hamiltonian 

h = h Q +n T.T-+ y"l. t + d Q Q 2 (69) 

in the basis of the 2s-ld shell states. 


The expression for the matrix elements of the 


Qq operator is as follows 


Cn’l'j'm’ |Q 0 |nljm> 


= <n'l ' |r 2 |nl> [ 3 2 3 ] L 3 2 3 ] 

10 1 
2 2 


m 0 m’ 

The Calcula t ion of Subsh e ll Occupation Numbers 


(70) 


The subshell occupation numbers for protons and 
neutrons were calculated as follows. The results of the 
HFB calculations are summarized in terms of the expansion 
coefficients C^ m and the amplitudes (U im »V ira ) appearing 
in eq (56). 

These values were used to first set up the 12X12 
matrix f n defined by eq (58). he next calculate the matrices 

(X p 

F ap( 0 ) * W ^ e - ) and defined by, respectively, the eqns 

(66), (65) and (64). These calculations were carried out 
for 20 values of the Gaussian quadrature points lying in 
the range (0,u/2)« Finally the subshell occupation 
numbers were calculated by using the eqs (60) and (63). 
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ASSUU'S, . .At© _ JJI SCUSSIG]] 

3~1 A Comparative Study of the PHF 3 and the Shell-iodel 
Approaches Employing the V'ildenthal Interaction 

Here we present the estimate for subshell 
occupation numbers for a number of 2s~ld shell nuclei in 
the frame work of the shell model as well as the PHF 6 
approaches, employing the 7/ildenthal interaction. The 
calculations are aimed at assessing the relative efficacy of 
the PHF 3 method vis-a-vis the calculation of occupation 
numbers. We have shown our results in Fig. 3 and Tables I 
and II. The results given here provide a first-ever explicit 
demonstration of the quantitative reliability of the PHF 3 
method. The PHF 3 results are throughout in excellent agree- 
ment with the shell model estimates; the maximum discrepancy 

32 

occurs in the case of the nucleus S where the i J RF3 method 
overestimates the shell model values by only about T.7^i. 

Vje have also given in Fig. 3 the experimental values 
for the subshell occupation numbers for neutrons resulting 
from various stripping and pickup reactions inoicated in 
Table I. The error bars associated with the experimental 
numbers displayed in Fig. 3 summarize the differences in 
the predictions resulting from various reactions. It is seen 
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that the Vildenthal interaction in conjunction with the 
shell model frame work provides a reasonably good overall 


description of the available experimental data. A noticeable 

P 8 

c’escrc-oancv however occurs in the case of the nucleus - Si 
where the Vilc’enthal interaction overestimates (underestimates) 
the subshell occupation numbers for the neutrons in the 
Id^/aCld^/p) orbits. This is not howev -r surprising in view 
of the fact that the structure of this nucleus involves 
maximum number of valance configurations for protons and 
neutrons in the 2s-ld shell. 

3-2 A Comparative Study of the PHF3 Predictions for Proton 
and Neutron Subshell Occupation Number Resulting from 
the Various Effective Interactions Operating in the 
2 s -lc 1 Shell. 

In Fig. 4 and Table II we have presented the 
subshell occupation numbers resulting x row various effective 
interactions described earlier. One notices wioely diff- 
erent predictions resulting from different interactions for 
nuclei heavier than 2 u\^g. The interaction labelled as K,Kd, 
K.+12F? and flip fail to reproduce the observed 
occupation numbers in the nuclei 2o iVg, 28 Si, 3 °Si» and o2 3. 

As far as the interactions I'.sSDI, ?b, and l are concerned 
they give nearly identical results agreeing reasonably well 
with the experiment throughout the shell for the (^ 5 / 2 ^ 

the ASDI interaction leads to extremely 


orbit 


vhereas 
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poor predictions for the (d 3 ^ 9 ) n orbit in the case of the 
PS 30 39 

nuclei 9 Si and S» the interaction P ! / leads to poor 

on 

results only in the case of S. 

In Table III we have compared the subshell occupa- 
tion numbers resulting from the PHF3 as well as the spectral 
method'*' for some N=Z nuclei in the 2S-ld shell. It is seen 
that the predictions resulting from the two methods are 
quite different. The differences are not entirely unanti- 
cipated since the spectral method is not expected to provide 
an accurate description of the properties associated with the 
ground states which are usually more than five widths away 
from the centroid of the energy distribution assumed to be 
Gaussian. 



TABLE I. Ground-state subshell occupation numbers 
from experimental data.^ 
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TABLE II. The subshell occupation numbers for protons 

and neutrons resulting from various effective 

interactions ara tabulated. The intrinsic 

quadrupole moments have been civon in units 
o 

of b , where b is the oscillator parameter. 





Protons Neutrons 


Interac- 

tions 

xr.2 v 

<o 0 > 


3 

2 

5 

2 

j= ! 

o 

o 

2 

5 

2 

K+SPE 

15.18 

0.53 

0.26 

1.30 

0.53 

0.26 

1.30 

K 3+SPE 

15.21 

0.60 

0.23 

1.30 

0.60 

0.23 

1.30 

K+12FP 

15.56 

0.67 

0.27 

1.15 

0.67 

0.27 

1.15 

RIP 

15. 30 

0.92 

0.15 

1.01 

0.92 

0.15 

1.01 

H5DI 

14.73 

0.58 

0.12 

1.38 

C. 58 

0.12 

1.38 

p:-; 

15.30 

0.60 

0.23 

1.26 

0.60 

0.23 

1.26 


15.36 

0.68 

0.20 

1.20 

0.68 

0.20 

1.20 

17 (oM) 


(0.51) (0.28) (1.21) 

(0.51) (0.28) (1.21) 

K+SPE 

18.94 

0 • 49 

0.56 

2.94 

0.49 

C. 56 

2.94 

K 3+SPE 

18.68 

0.50 

0.43 

3.07 

0.5 C 

0.43 

3.07 

K+12FP 

19.01 

0.42 

0.66 

2.90 

0.42 

0.66 

2.90 

RIP 

18.98 

0.68 

0.37 

2.99 

0.68 

0.37 

2.99 

jViSDI 

17.69 

0.50 

0.18 

3.36 

0.50 

0.18 

3.36 

pw 

17.61 

0.26 

0.28 

3.49 

0.26 

0.28 

3.49 

0 

18.03 

0.33 

0.33 

3.38 

0.33 

0.33 

3.38 

L(Siu) 


(0.45) (0.56) (2.99) 

(0.45) (0.56) (2.99) 

K+SPE 

20.98 

0.45 

0.38 

3.21 

0.60 

1.77 

3.63 

K 3+SPE 

-18.28 

0.98 

0.37 

2.52 

1.47 

0.78 

3.75 

K+12FP 

21.03 

0.38 

0.47 

3.18 

0.43 

1.96 

3.61 

RIP 

-18.31 

1.48 

0.25 

2.21 

1.68 

0.85 

3.47 

Jr.SDI 

-15.33 

0.76 

0.12 

2.98 

0.95 

0.20 

4.85 

Pi/ 

-15.71 

0.34 

C.27 

3.26 

0.64 

0.39 

4.97 

II 

-15.99 

0.43 

0.27 

3.12 

0.73 

0.39 

4.88 


W(,SM) (0.35) (0.45) (3.20) (0.56) (0.62) (4.82) 


Contd 
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Contd... TABLE II 


Nucleus Interac- 
tions 


< Q o> 


Protons 

1= k 1 
J 2 2 


5 

2 


Neutrons 
1 3 5 

1- 2 2 2 


28 


Si 


'Si 



K+Sl-S 

-22.77 

K3+SPE 

-22.24 

N+121-P 

-23.01 

RIP 

-22.37 

MSDI 

-18.85 

PW 

-19.82 


-20.28 

i/(sm) 

n+sph 

-20.01 

KB+SPE 

-18.65 

K+12FP 

-20.46 

RIP 

-19.78 

MSDI 

6.96 

PYJ 

9.76 

v/ 

-15.69 

V' ( Sil) 

K+SPE 

-17.52 

K3+SPE 

-13.62 

K+12FP 

15.03 

RIP 

-16.91 

MSDI 

0.43 

py; 

0.08 

w 

13.21 

W(SN) 


1.33 

1.12 

3.55 

1.45 

0.97 

3.58 

1.26 

1.25 

3.49 

1.65 

1.09 

3.25 

1.12 

0.37 

4.51 

0.90 

0.61 

4.49 

0.99 

0.65 

4,36 


(0.70) (0.68) (4.62) 


1.29 

0.36 

3.85 

1.45 

0.66 

3.88 

1.13 

1.01 

3.86 

1.71 

0.93 

3.35 

1.05 

0.11 

5.02 

0.10 

0.06 

5.83 

0.59 

0.24 

5.17 


(0.68) (0.58) (4.74) 


1.36 

1.70 

4.98 

1.78 

1.00 

5.27 

0.23 

2.14 

5.63 

1.77 

1.61 

4.69 

1.98 

0 . 04 

5.94 

1.98 

0.04 

5.94 

1.14 

0.96 

5. 86 


(1.41) (1.16) (5.42) 


1.33 

1.12 

3.55 

1.45 

0.97 

3.58 

1.26 

1.25 

3.49 

1.65 

1.09 

3.25 

1.12 

0.37 

4.51 

0.90 

0 . 61 

4.49 

0.99 

0.65 

4.36 

(0.70) (0.68) ( 

[4.62) 

1.41 

1.87 

4.78 

1.64 

1.55 

4. 89 

1.26 

2.15 

4,76 

1.67 

1.82 

4.55 

1.95 

0.07 

5.94 

0.58 

1.45 

5.93 

1.22 

1.29 

5.48 


(1.41) (1*44) (5.15) 


1.36 

1.70 

4.98 

1.78 

1.00 

5.27 

0.23 

2.14 

5.63 

1.77 

1.61 

4.69 

1.98 

0.04 

5.94 

1.98 

0.04 

5.94 

1.14 

0.96 

5.86 


(1.41) (1.16) (5.42) 


Contd. . . . 





Contd. . . . 

T AISLE 

II 
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Protons 


Neutrons 

Nucleus 

Interac- 

tions 

• <> 

> 1 

O 

2 

5 

2 

. 1 

J= 5 

O 

O 

2 

5 

2 

34 s 

K+SPE 

-14.75 

1.32 

1.40 

5.23 

1.57 

2.61 

5.75 


KB+-SPE 

-10.01 

1.88 

0.42 

5.67 

1.94 

2.12 

5.86 


K+12FP 

-14.54 

0.47 

1.96 

5.52 

0.68 

3.48 

5.77 


RIP 

-12*73 

1.88 

1.15 

5.04 

1.89 

2.14 

5.89 


iviSDI 

1.21 

1.98 

0.04 

5.94 

1.97 

2.08 

5.94 


PW 

-5.85 

1.97 

0.05 

5.93 

1.98 

2.01 

5.94 


v i 

\D 

00 

* 

1 — 1 
rH 

1 

1.46 

0.70 

5.77 

1 . 66 

2.38 

5.88 


'V ( Sit) 


(1.66) (0.73) (5.61) 

(1.76) 

(2.48) (5.76) 

36 Ar 

K+SPE 

-14.37 

1.60 

2.64 

5.68 

1,60 

2.64 

5.68 


K3+SPu 

-11.46 

1.9C 

2.19 

5.82 

1.90 

2.19 

5.82 


11+1 2 FP 

-15.34 

0.93 

3.31 

5. 68 

0.93 

3.31 

5.68 


PJ? 

-10.16 

1.93 

2.08 

5.91 

1.93 

2.08 

5.91 


NSDI 

-11.65 

1.84 

2.18 

5.90 

1,84 

2.18 

5.90 


PW 

-9.52 

1.94 

2.05 

5.93 

1.94 

2.05 

5.93 


w 

-13.52 

1.62 

2.4 5 

5.84 

1.62 

2.45 

5,84 


'k ( SA-.) 


(1.78) (2.68 

) (5.54) 

(1.78) (2.68) (5.54) 

The 

numbers 

given in 

brackets denote the 

exact 

shell 

model 

results obtained wit] 

n the 

Nildenthal (W) 

interaction. 




46 


TABLE III. Comparison of the subshell occupation 

numbers resulting from the PHF3 as well as 
the spectral method for some M=2 nuclei 
in 2s -Id shell. 


Spectral Method K'lFB Method 


Nucleus 

Interaction 

i=§ 

3 

2 

5 

2 

i4 

3 

2 

5 

2 

~°Ne 

P T 7 

0.5 

0.1 

3.4 

1.2 

0.5 

2.5 


K+12FP 

0.6 

0.2 

3.1 

1.3 

0. 5 

2.3 


HI? 

1.7 

0.0 

2.2 

1.8 

0.3 

2.0 


i -;sdi 

0.5 

0.0 

3; 5 

1.2 

0.2 

2.8 

24 „ 

p V ; 

0.3 

0.5 

6.7 

0.5 

0 . 6 

7.0 


K+12FP 

0.9 

1.2 

5.9 

0.8 

1.3 

5.8 


RIP 

3.2 

0.4 

4.5 

1.4 

0.7 

6.0 


/‘•:sdi 

i;5 

0.1 

6.4 

1.0 

0.4 

6.7 

28 Si 

Plv 

1.4 

1.1 

9.5 

1.8 

1.2 

9.0 


K+12FP 

1.4 

2.3 

8.3 

2.5 

2.5 

7.0 


RIP 

3» 6 

1.2 

7.2 

3.3 

2.2 

6.5 


7 IS LI 

2; 3 

0.5 

9.2 

2.2 

0.7 

9.0 

32- 

P'7 

2.7 

2.0 

11.3 

4.0 

0,1 

11.9 


R+12r r 

1 O 

3.8 

10.4 

u. 5 

4.3 

11.3 


RIP 

3.9 

2.4 

9.7 

3.5 

3.2 

9.4 


M-SDI 

3.3 

1.4 

11.2 

3.9 

0.1 

11.9 

O < 

° 3 Ar 

PR 

3.8 

4.2 

12.0 

3.9 

4.1 

11.9 


K+12FF 

2.3 

6.0 

11.7 

1.9 

6.6 

11.4 


RIP 

4.0 

4.3 

11.7 

3.9 

4.2 

11.8 


MSB I 

3.9 

4.1 

12.0 

3.7 

4.4 

11.8 




FIG-1 T= 1 MATRIX ELEMENTS OF VARIOUS EFFECTIVE INTERACTIONS 







FIG. 2 T= 0 MATRIX ELEMENTS 




SHELL MODEL 
PHFB 



FIG. 3 COMPARISON OF THE RESULTS OBTAINED IN THE SHELL MODEL AND THE PHFB 
FRAMEWORK 
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FIG. 4 SUBSHELL OCCUPATION NUMBERS FOR NEUTRONS VIA THE PHFB METHOD 
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